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WITH A SLENDER FUSELAGE AT HIGH SUPERSONIC SPEEDS 
By Thomas N. Canning ajid Billy Pat Denardo 

SUMMARY 


Tests were conducted In the Ames supersonic free-flight wind tunnel 
to evaluate the usefulness of available theory for the calculation of 
lift -curve slope and center of pressure of low -aspect-ratio, triangular, 
and rectangular cruciform wings in combination with a long slender body. 
Tests were made in the range of Mach numbers between 1.3 sind 6.2 and at 
Reynolds numbers from 2.8 million to l6.0 million based on model length. 
The results show that theoretical calculations give the lift-curve slope 
near zero lift within the ejqjerimental scatter and the center of pressure 
within about 3 percent of the body length. The variation of minim\im drag 
coefficient with Mach number is presented for both models. Estimates 
based on experimental values of lift-curve slope and drag indicate that 
the triaingular wing-body combination is the more efficient and may develop 
a. lift-drag ratio of 4.5 at a Mach number of 6.0 for power-on operation, 
neglecting body-base drag. 


INTRODUCTION 


Many tests of wings and wing-body combinations have been made and 
the resxilts compsired with linear and more exact theories in the Mach num- 
ber range below 3* In the Mach number range above 3, on the other hand, 
there are very few congiarisons available (references 1 a nd 2) . The pres- 
ent investigation is Intended to provide some data for low -aspect -ratio, 
triangular, and rectangular cruciform configurations in this speed rauige 
at Reynolds niimbers corresponding to a 10-foot-long missile flying between 
80,000 and 100,000 feet altitude. The usefulness of available theory is 
investigated by con 5 >arlng the calculated and experimentally determined 
aerodynamic coefficients in the range of Mach numbers from 1.3 to 6.2. 

The rsuige of Reynolds numbers based on body length was 2.8 million at 
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M = 1.3 to l6.0 million at M = 6.0. The effect of Reynolds number on 
the aerodynamic coefficients was investigated at M = 6.0. 


syMBOIS 
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horizontal projection of the instantaneous acceleration of 
the model center of gravity normal to the tunnel center 
line, feet per second squared 

(*) 

drag coefficient at zero lift 

(4) 

pitching-moment coefficient about model center of gravity 


drag coefficient 


lift coefficient 


I pitching moment A 

V y 


/dC- 


lift-curve slope (^“5a) * degree 


pitching -moment -curve slope degree 

drag, pounds 
Naperian base 

constants defining variation of a with time 


frequency of the pitching motion, cycles per second 

principal moment of inertia about lateral axis through model 
center of gravity, foot-pound seconds squared 


ratio of two-dimensional lift-curve slope of airfoil calcu- 
lated with shock-expansion theory to that calculated with 
linear theory 


damping constant, per second 
model length, feet 


L lift, pounds 

m model mass, slugs 


MCA RM A52C24 


3 


M test Mach maniber 

cLq free-stream dynamic pressure, poimds per foot squared 

Q double integral with respect to time of the angle of attack 

function, seconds squared 

R Reynolds number based on model length and free -stream 

conditions 

S maximum cross-sectional area of model body, feet squared 

t time, seconds 

V horizontal projection of the instantaneous velocity of model 

center of gravity normal to tunnel center line, feet per 
second 


^C.g.-'^c.p. distance from model nose to center of gravity and center of 
pressiire, respectively, feet 

y horizontal distance from t unn el center line, feet 

a. angle of attack of model relative to local flight path, 

degrees 

w 2nf 


Subscripts 


1 , 2 , 3 , 4 refer to times, positions, and velocities at the instants of 

exposure of the shadowgraphs in stations 1, 2, 3, and 4, 
respectively 

12,14, refer to intervals between stations 1 and 2, 1 and 4 , 

23, etc. 2 and 3, etc. 


EXPERIMENTAL TECHNIQUES AND MODELS 
Facility 


The experiments were performed in the supersonic free-fll^t wind 
tunnel of the Ames Aeronautical Laboratory. This facility is a short 
ballistic range within a variable pressure, supersonic, blowdown wind 
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tiinnel. Details of timnel design and operation are given in reference 3. 
The features of the facility pertinent to the present Investigation are 
mentioned below. 

The model is fired from a smooth-bore gun located in the wind- 
tunnel diffuser. As it passes upstream through the 15-foot-long test 
section 7 shadowgraph plctiires, 4 in the horizontal plane and 3 in the 
vertical plane, are made. A chronograph records the instants of exposure 
of the shadowgraphs, completing a time history of the model position. 

Two typical shadowgraphs are shown in figure 1. 

Methods of Measuring Aerodynamic Coefficients 


General description.- The model is designed to execute between one - 
half <~>np and one -quarter pitching oscillations in the test section 
is disturbed at launching so as to oscillate with an amplitude of 
about 5° in the horizontal plane. Small oscillations in the vertical 
plane result from accidental disturbances. The complete motion in 
3 dimensions cannot be studied because of inadequate data in the vertical 
plane. Instead, the projection of the motion in the horizontal plane is 
used, assijmlng that the interaction between pitch and sideslip is negli- 
gible. 

The aerodynamic coefficients are calculated from the record of model 
motion as a function of time. For an example of such data, see fi^e 2. 
Two methods have been used for determining lift-curve slope, pitching- 
moment-curve slope, and center-of -pressure position, 0^^, Cm^, Xc.p., 

from the behavior of models oscillating in pitch. One method consists of 
measuring the flight-path curvature due to lift and the other consists of 
measuring the change in pltchlng-moment-c\irve slope resulting from large 
center-of -gravity shifts . 

Lateral -movement method.- This method uses flight -path ctirvature as 
a measure of lift. Using Newton’s second law of motion, the instan- 
taneous acceleration of the center of gravity of the model normal to the 
tunnel axis may be written as follows: 


a 


= ^ = 
dt^ 


CLq, O' <loS 
m 


( 1 ) 


This equation assumes a linear lift curve and neglects the contributions 
of the drag force, and the lift force due to pitching and plunging. 

These contributions are usually, but not always, negligible so care was 
exercised to see that the omissions were permissible. The direction of 
the lift force is, by definition, perpendicular to the flight path, but 
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the maxinrum observed inclination of the flight path relative to the 
tunnel axis vas 0.5^. The error in the experimentally determined lift 
introduced here is given by the cosine of 0.5^ and is therefore small* 


The time variation of angle of attack in equation (l) can be obtained 
using the angle -of -attack measurements from the shadowgraph pictures and 
two assumptions which define the form of the motion. The assumptions are 
that the restoring moment is directly proportional to the angle of attack 
and the damping moment is directly proportional to the pitching and 
plunging rates. This leads to a variation of angle of attack given by 


a 


-kt 

e 


(E cos cot - F sin cot) 


( 2 ) 


Equation (2) is fitted to the observed variation of a with respect to 
time by a least-squares procedure described in reference k. In this way, 
the four unknowns, k, E, and F, are evaluated. The complete set of 
4 shadowgraphs is required to determine the sine wave so the method of 
least squares is not strictly required but is used as a convenient and 
systematic method. 

Combining equations (l) and (2) gives the instantaneous lateral 
acceleration as a function of time: 


^ e“^^ (E cos cot - F sin cjt) (3) 

dt dt m 


Integrating with respect to t gives the following equation for lateral 
velocity as a function of time: 


Y 


dy r* 

37 ' * -^Jt, 


(E cos u)t - F sin o)t)dt 


(M 


in which Vi is the lateral velocity at the first station. A second 
integration gives the equation of lateral displacement. The integrals 
are evaluated between limits corresponding to the times and lateral 
positions in two shadowgraph stations. 


72 = 7i + Vi(ts 



e“^^(E cos o)t - F sin ut)dt dt 

(5a) 


The double integral on the right can be evaluated since both the integrand 
and the limits are known. For brevity, it will be designated 


^2 = yi + 




(5to) 


m 
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In this equation, the unknowns axe and Vi. The constant 

be eliminated by using data from the third station: 


^3 = * ''i(*3 - + 




m 


Q 


13 


cam 


(5c) 


Solving equations (5b) and (5c) simultaneously yields the following 
expression for C]^: 

ni yjg ' yi3(^ia/^i3) 

^12 “ ^13(^12/^13) 


The least-squares fit to the variation of angle of attack with time 
establishes the pitching frequency f. This makes it possible to obtain 
the pitching -moment -curve slope about the center of gravity from the 
relation 


f 


1 

" 2rt V I 


or 



f^ I 
SqoZ 


(7) 


The center of pressure can be obtained from the lift and pitching-moment 
results using the relation: 


( 





(8) 


Pitching -moment method ,- The relation between lift-curve slope^ 
center of pressure, and pitching -moment -curve slope given In equation (8) 
suggests that tests be made of slmllax models with widely separate centers 
of gravity. Solving equation (8) simultaneously for two such tests deter- 
mines the unknowns Ci^ and Xc.p. since for each test Is given by 

equation (T)- This method has also been used In this Investigation. 


The lateral -movement method Is Inherently more accurate than the 
pltchlng-moment method for this facility and yields much smaller scatter 
In Both methods give reliable results for center of pressure, but 


the lateral -movement -method results contain less scatter here also. 
Values of obtained with both methods sire presented subsequently, 

but the values of Xc.p./^ from only the lateral -movement method are 


presented. 


Drag . - From the time history of the model position, deceleration and 
hence drag were determined. These calculations are based on the 
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assvtmption that Cj) is constant. Details of this part of the data 
analysis are included in reference 3* Since, in these tests, angles of 
attack larger than 8.0° were sometimes observed, the varying drag due to 
lift was important. The effect of this varying drag, although not 
treated exactly, was accounted for approximately by subtracting from the 
indicated average value of Cj) an approximate value of drag due to lift 
given by the mean value during test of (l/57*3) X a^. 


Models 


The models were wing-body combinations and were constructed of 
various metals chosen to give the desired wei^ts and center-of -gravity 
positions. The bodies were cone cylinders of fineness ratio 15 with a 
nose length of 6 body diameters. Wings of triangvilax and recteingular 
plein form were used. The aspect ratio of the exposed-wing panels joined 
together was 0.64. The profiles were of constant thickness except for 
bevelled leading edges; the constant thickness resiilted in a constant 
airfoil section of 0.051 thickness ratio for the rectangular wing and a 
thickness ratio varying from 0.026 at the root to 0 .l 60 at the tip on the 
triangular wing. This design compromise was accepted to reduce cost of 
the models, and it is believed that the lift and center-of-pressure 
results would not have differed significantly if the same airfoil section 
had been used on both models. Detailed sketches of the models are shown 
in part (a) of figure 3 and photographs are presented in figure 4. 

The models were fired from a smooth-bore 20-mm. gun and were held 
in plastic sabots designed to give support and to inpart a pitching 
distiirbance to the model at separation. A sketch of a typical sabot is 
given in figure 3(l3)^ and a photograph is presented in figure 4. The 
pitching motion resulting from this initial distirrbance was recorded as 
described above to give the data reqxiired. As may be seen in figure 3(i>) 
and figure 4, the models left the gun with an initial angle of attack, 
and the plane of this angle was parallel to the plane of the 4-shadowgraph 
group. 


THEORETICAL CALCDLATIOrJS 


Theoretical values of lift-curve slope and center of pressure of 
each configuration were calculated at several points within the Mach 
number range. An approximate theory similar to that of Nielsen 
(references 5 and 6) was iised for both configurations. In the reference 
papers^ the approach is to combine the calciilated lift of the various parts 
of the model alone and the lift induced on each part by the others. The 
load distribution on the body alone is calculated using slender-body 
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theory. The lift distribution on the two exposed wing panels joined 
together is estimated using Stewart’s paper (reference 7) for the 
triangiaar wing and the paper by Lagerstrom (reference 8) for the low- 
aspect-ratio rectangular wing. The lift induced on the wings by the body 
upwash field is estimated using the slender wing-body theory of Spreiter 
(reference 9). The point of application of this lift component is taken 
as the center of pressure of the two wing panels joined together. The 
lift induced by the wing on the body and its point of application are 
estimated by reducing the problem to the planar case in a manner similar 
to that developed by Morikawa (reference 10) . Superposition of these 
components of lift and pitching moment gives the desired estimates of 
and Xc.p./I* 

In the present investigation, a modified theory was used in which 
foiir departures from the method as described were made, as sioggested in 
reference 5, in order to make the theoretical flow correspond more 
closely with the real flow . The lift on the nose of the body was calcu- 
lated using Kopal's tables (reference 11), and the lift carry-over from 
the nose onto the cylindrical portion of the body was evaluated by means 
of Tsien's paper (reference 12), using the exact rather than the approxi- 
mate relation between local velocity and pressure. In calculating the 
lift carry-over onto the body from the wing, the fact that these models 
have no body aft of the wing was accounted for approximately . In addi- 
tion, in the case of the rectangular wing, the sum of lift on the 
rectangular wing and the lift on the body due to the presence of this 
■ying was increased by accounting for the fact that the flow consisted of 
a shock-wave expansion flow instead of the Mach wave type of flow assumed 
in reference 6. This was done by mixltiplying this part of the lift by 
the ratio 

Ct S.E. 

K = -is 

Ct L.T. 

in which Cj^S.E. is the lift -curve slope of the actual profile in two- 
dimensional flow calculated using shock-expansion theory, and Cj^L.T. is 

the same quantity calculated using linearized theory. This procedure was 
used in reference 2. The theoretical contribution of each part of the 
configurations to the lift is presented in figure 5* fh.® iiet effect on 
the estimated lift-curve slope of the rectangular-wing model introduced 
by the four modifications is large at high Mach numbers. 


RESULTS AND DISCUSSION 

The experimental variation of Cj^ with Mach number for the two 
models is given in figures 6 and 7 for both types of analysis of the same 
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test firings. Comparison of the two sets of results indicates a fair 
consistency between the two methods except that the lateral -movement 
method yields somewhat lower values of at the higher Mach nijmbers 

and yields less scatter at all Mach numbers. For purposes of discussion, 
the results of only the lateral -movement analysis will be considered 
because of the apparent advantage in accuracy. The small n umb er of data 
points on these plots resiilts from the requirement that displacements 
from a straight flight path of O .15 inch or greater must be observed in 
order to obtain reliable answers. One possible explanation of the 
remaining scatter in these results is described in the appendix. The 
values of predicted using the modified theory were within the 

apparent band of experimental scatter of ±9 percent from the median value 
for both models at all Mach numbers at which tests were made. This indi- 
cates that existing theoretical results are of value in this Mach number 
range for simple configurations. 

The experimentally determined centers of pressure, obtained by the 
lateral -movement method, were usually forward of the predicted positions 
as shown in figure 8 . The maximum discrepancies between experiment and 
theory were observed on one model of each type at M = 4.5; in each case 
the discrepancy was about 5 percent of the body length. The forward 
movement of the center of pressure with increasing speed was well 
predicted. 

The tests show that the rectangular wing exhibits slightly superior 
effectiveness as a stabilizer for the test body both because its center 
of pressiire is farther aft and because it has a higher lift-curve slope. 
This superiority is gained at the expense of a large drag increase over 
that of the triangular wing as is shown in figure 9* This figure 
presents the variation with Mach niomber of the zero-lift drag, . , 

of both models. ^ 

Although direct measurement of lift and drag at constant angle of 
attack was not possible in this investigation, a comparison of lift-drag 
ratios was made based on the assiamptions that 

Cl = Ci^a 


and 


% = c 


D, 


I + ^ Ct a‘ 

min 57*3 


neglecting the lift and drag of whatever trimming device might be chosen. 
The angles of attack for maximiam lift-drag ratios at M = 6.0 were 
7 . 5 ° or less. Since the amplitude of oscillation of the test models was 
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usually about 5*0^^ it is felt that extrapolation to T»5 should not 
introduce important errors in the comparisons. When the value of 

was taken to be the total measured value at zero lift, the attainable 
lift-drag ratios of the two models at M = 6.0 are ^v.O and 3-^ for the 
triangular and rectangular wing models, respectively, at corresponding 
lift coefficients of 1.^ and 1.7. 

To estimate the performance which might be realized in power-on 
operation, the drag was reduced by subtracting the estimated body-base 
drag (base pressure assumed equal to 0.3 of the free-stream static 
pressure) . This assumption is conservative as judged from the data of 
reference 13 although the Mach nvimber range of the reference test was 
somewhat lower. The triangular-wing model again appeared more efficient 
having, at M = 6.0, an (L/D)niax of at Cl = 1-2, as compared to an 

(l/d) of 3.6 at Cl = 1*6 for the rectangular -wing model. These 
values^f attainable lift-drag ratios are considered encouraging in view 
of the fact that extreme aerodynamic efficiency was not the basis for 
design. The comparative efficiencies of the two models might be altered 
somewhat by using the same airfoil section on both models, but the changes 
from the present results would probably be small. 

The results of tests indicated no measiirable effect of Reynolds 
number on the lift-curve slope, center of pressure, or drag of the 
models at M = 6.0 in the range of Reynolds number from 5*3 to l6.0 
million. This corresponds to a Reynolds number range based on mean 
exposed chord of 1.2 to 3*6 million. 

CONCLUDING REMARKS 


The theory of Nielsen (references 5 smd 6) , when altered to account 
for the effects of airfoil section, lift carry-over from the nose onto 
the cylindrical portion of the body, and the absence of an afterbody, 
predicts values of Ci^ within the experimental scatter and the mean 
center-of -pressure positions within about 3 percent of the body length 
for both models at Mach nximbers between 1.5 and 6.0. 

Although the rectangular wing is slightly more effective than the 
trieingular wing as a stabilizer for the test body, estimates of the lift- 
ing efficiency of the two models indicate that the triangular wing is 
superior as a lifting surface. Estimates of the lift-drag ratios for 
power-on operation with no allowance for the effects of a trimming device 
indicate that values as high as U.5 may be attained at a Mach niimber of 6.0 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX 

APPARENT SCATTER IN Cj_^ 


The method of data analysis described in the body of the report is 
based on the assumption of linear pitching -moment characteristics. 
Experimental results for many similar configurations indicate seriously 
nonlinear characteristics at large angles of attack. Since the present 
investigation was restricted to small angles, nonlinearities probably 
did not contribute important errors in the results^ but may have intro- 
duced some scatter. Such pitching -moment characteristics are frequently 
well represented by the expression 

^m “ ^m(x^ constant 

Assuming this in the present case, the oscillatory motion resulting from 
an initial disturbance is complicated, but it may be seen that since the 
restoring moment is disproportionately large at high angles of attack 
the peaks will be sharper than in the case of a sine wave of the same 
frequency and amplitude. 

Assume that the actual oscillation is given in the sketch by the 
solid line 
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Consider first that the positions of the shadowgraph images correspond 
to the circled points. If a sine wave is passed through these points, 
indicated by the dashed curve, it is seen that the area under each loop 
of the sine wave is significantly larger than the area under the actual 
curve. Since the lift-cxxrve slope is obtained by comparison of the 
lateral movement of the model to the integration of this curve twice 
with respect to time, as described in the body of the report, it may be 
seen that the values of Q in equation (6) will be too large and the 
analysis will underestimate the average lift-curve slope. On the other 
hand, if the shadowgraphs correspond to the ^ points, the lift-curve 
slope will be overestimated. 

The results of an examination of all test runs for conditions at 
which duplicate data was obtained support the above contention. Two 
tests which fitted the pattern for underestimation, including each type 
of model, gave decidedly low answers. Similarly, three tests which 
fitted the pattern for overestimation gave decidedly high answers. Nine 
tests, which fitted neither pattern, either agreed well with one another 
or gave results which fell inside the scatter of other runs. 

Three tests did not agree with the correlation. Two of the three 
did not fit either pattern and gave high answers. The amplitude of 
pitching oscillation in one of these cases was small and the accuracy of 
measuring y may not have been sufficient; the pitching amplitude in the 
other case was greater than 10° and may have gone into the seriously 
nonlinear lift-curve range. The third test fitted the pattern for under- 
estimation and gave a high answer. The pitching amplitude in this last 
case was small enough to indicate that accuracy of measuring y may 
again have been marginal. 

It is believed that the errors introduced by the nonlinearities 
discussed above do not seriously impair the value of the experimental 
results. The median point of the scatter probably represents quite well 
the average value of lift -curve slope over the experimental range of 
angle of attack. 
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(a) Triangular-wing model, a = -3.0°. 



Distance from tunnel center tine, feet Angle of attack, degrees 
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(a) Variation of angle of attack with time. 
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(b) Variation of distance from tunnel center line to mode! center 

of gravity with time. 

Figure 2. -Time history of mode! motion. 
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Section A- A 




(a) Models. 
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Figure 3,-Test configurations. 
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A-16856.1 

(a) Models and disassembled sabot. 



(b) Partially assembled test round. 

Figure k.- Photographs of test configurations. 
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(a) Latera/- movement method. 

Figure 6.-Voriotlon of with Moch number for trionguiar-wing mode/. 
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fdJ Pitching- moment method. 
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(a) Lateral- movement method. 


Figure 


7.-Variation of 


C, with Mach number for rectanguiar-wing model. 
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(b) Pitching- moment method. 


Figure 7. - Concluded. 


ro 


NACA RM A52C2U 


26 


NACA RM A52C24 




— djnssdjd yo JdiUdQ 

dSOU UJOJJ. 


(a) Triangular wing. 

Figure 8.-Variation of center-of -pressure location with Mach number. 
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(b) Rectangular wing. 
Figure 8. - Concluded. 
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Figure 9. -Variation of drag coefficient with Mach number. 


